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SECTION  I 
INTRODUCTION 

The  purpose  of  this  paper  is  to  characterize  scintillated  UHF  signal  struc- 
tures in  expected  propagation  environments  for  either  ambient  or  nuclear  envi- 
ronments, The  characterization  of  UHF  scintillation  effects  is  necessary  for 
the  evaluation  of  existing  and  proposed  UHF  satellite  communications  systems. 
UHF  systems  are  likely  to  be  very  important  for  many  years  to  come  for  several 
reasons.  First,  there  already  exists  a great  deal  of  UHF  operational  experi- 
ence from  which  to  draw  in  developing  new  systems.  Second,  UHF  hardware  tech- 
nology is  well  developed,  lowering  system  costs.  Also,  propagation  through 
disturbed  environments  requires  high  rf  power  levels  which  favors  lower  fre- 
quencies with  their  higher  dc  to  rf  power  conversion  efficiency.  Finally,  UHF 
systems  appear  to  be  adequate  for  many  important  satellite  communication  appli- 
cations despite  their  higher  sensitivity  to  various  propagation  effects. 

This  paper  consists  of  five  sections.  The  next  section  defines  the  propa- 
gation environments  used  to  examine  the  basic  characteristics  of  UHF  propaga- 
tion. The  third  section  discusses  the  propagation  methods  that  are  used.  The 
fourth  section  presents  the  results  of  the  propagation  calculations  and  the 
implications  regarding  UHF  systems.  The  final  section  will  summarize  the  basic 
conclusions  of  this  work. 
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SECTION  II 

PROPAGATION  ENVIRONMENTS 


The  propagation  environment  is  a layer  containing  electron  density  fluctua- 
tions lying  perpendicular  to  the  propagation  line  of  sight  (LOS).  Figure  1 por- 
trays the  geometry  and  the  relevant  dimensions.  The  values  for  and  z^  were 
chosen  to  be  a reasonable  compromise  between  nuclear  and  ambient  environments. 
The  value  of  z^  may  be  realizable  in  nuclear  environments  and  can  be  approached 
at  low  elevation  angles  in  ambient  environments. 

The  fluctuation  statistics  in  the  layer  are  isotropic  in  two  dimensions, 
one  of  which  is  parallel  to  the  LOS.  No  variation  in  electron  density  is  as- 
sumed in  the  third  dimension.  The  fluctuation  statistics  are  uniform  over  the 
layer  and  are  described,  in  part,  by  one  of  four  correlation  functions.  The 
remainder  of  the  statistical  description  is  supplied  by  either  assuming 
gaussian  fluctuations  or  by  using  some  explicit  representation  for  the  fluctua- 
tions. Both  the  gaussian  assumption  and  explicit  representations  will  be  used 
to  demonstrate  that  for  sufficiently  large  z^  only  the  correlation  functions 
or,  equivalently,  the  second  order-statistics  are  necessary  for  calculating 
propagation  effects. 

The  electron  density  fluctuation  correlation  functions  are  characterized 
by  the  mean  square  electron  density  fluctuation,  An2,  an  outer  scale  size,  L, 
and  possibly  an  inner  scale  size,  i.  Except  where  otherwise  noted  the  values 
of  the  outer  and  inner  scale  sizes  will  be  1.25  km  and  30  meters  respectively. 
The  assumed  value  for  L is  the  smallest  measured  length  associated  with  high 
relative  fluctuation  power  (Anz/n^  = o.8)  known  to  the  authors.  The  value  for 
the  inner  scale  was  chosen  because  measurements  with  about  30  meter  resolution 
in  scale  size  in  the  ambient  ionosphere  did  not  reveal  an  inner  scale. 

The  first  correlation  function  used  to  represent  the  electron  density 
fluctuations  is  exponential  in  nature.  A simple  candidate  form  is 


R|^-2(e,p) 

p 

R|^-2(e,p) 

P 


n(z+e,x+p)  n(z,x) 


Anz 


(la) 

(lb) 
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This  case  will  be  labeled  as  because  the  power  spectral  density  correspond- 
ing to  equation  (lb)  is  proportional  to  K“2  for  > 1/L.  Power  spectral  densi- 
ties with  this  dependence  are  characteristic  of  ambient  environments  as  seen  in 
data  from  in  situ  measurements  (refs.  1 through  4)  and  calculations  (refs.  5 and 
6).  Equation  (1)  is  usually  used  when  the  fluctuations  are  taken  as  gaussian 
distributed  and  i « L.  When  an  explicit  representation  of  n(z,x)  is  required, 
the  fom  in  equation  (lb)  is  slightly  modified  but  retains  the  exponential 
nature. 

To  construct  an  explicit  representation  let  us  consider  a group  of  randomly 
located  gaussian  rods  of  electron  density  in  some  area  A where  A is  in  a plane 
containing  the  LOS  and  other  fluctuation  coordinate.  Let  us  further  assume 
that  there  are  on  the  averabe  fi/A  rods  distributed  per  unit  area.  The  two- 
si  onal  electron  density  is  then 


where  a^.  and  N are  distributed  according  to 


P;,(a.)  = t/v^a?,  a.  > i/^ 

a I II 


p (N)  = rL® 


4An^  A/  NllLir3/2 


(l-erf  {|- 


1/2  ^-a2/4L2 


where 


erf(x)  - ~ I dt 

^ / - 


Note  that  now  both  an  inner  and  outer  scale  size  is  required. 
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After  some  manipulation,  it  can  be  shown  that  the  one-dimensional  correla- 
tion function  and  power  spectral  density  are,  respectively. 


R^-2(e) 

p 


PSD(Kp)  = 


2An2e 

L 


-£2/4L2 


-K2£2/4n 

_e ^ 

1/L2  + k2 

p -■ 


(7) 


The  mean  electron  density  is 


7r3/2N  n (O)ilL  r 
nTF,77  = [l-erf 


(8) 


when  «,  « L,  equations  (6)  and  (7)  demonstrate  the  required  statistical  proper- 
ties. It  is  also  interesting  to  note  that  neither  the  power  spectral  density 
or  the  correlation  function  depend  on  N/A.  This  means  that  F/A  can  be  varied 
independent  of  the  second  order  statistics.  Intuitively,  as  U/k  increases,  the 
number  of  rods  intersected  by  a straight  line  ray  must  increase  or,  conversely, 
the  mean  free  path  between  intersectionr  must  decrease.  Let  us  define  a mean 
free  path  (MFP)  as 


[2a][^a/L][e-""^^‘-"]  P^la)^^ 


(9) 


This  definition  is  analagous  to  conventional  scattering  theory  where  the  inte- 
gral represents  the  average  scattering  cross  section  weighted  for  the  differing 
size  and  magnitude  of  the  gaussian  rods.  A MFP  will  be  used  as  the  measure  of 
thickness  of  scattering  layers  to  determine  when  only  the  correlation  function 
is  sufficient  to  do  the  propagation.  Evaluating  the  integral  in  equation  (g), 
we  have 
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MFP  E 


2Wl^ 


1 - erf 


-1 


(10) 


The  explicit  representation  just  described  may  or  may  not  resemble  the  actual 
environment  for  higher  than  second  order  statistics.  The  representation  was  de- 
veloped only  to  study  the  dependence  of  the  propagation  on  the  higher  order 
statistics  by  varying  the  MFP  layer  thickness  or  by  comparing  with  the  gaussion 
fluctuation  assumption. 

_2 

The  second  correlation  function  to  be  considered  is  the  inverse  of  the  Kp 
case,  that  is,  the  correlation  function  is  now  power  law.  This  will  be  re- 
ferred to  as  the  intermediate  case. 


*^INT  " An2/(1  + (e2+p2)/2L2)  3/2 


(11) 


The  gaussian  assumption  on  the  fluctuations  will  always  be  used  in  this  case. 
This  case  has  been  postulated  for  some  types  of  disturbed  environments. 

The  third  correlation  function  used  is  a gaussian  function. 


Rg(e,p)  = An2 


(12) 


This  function  is  not  necessarily  thought  to  be  characteristic  of  any  type  of 
environment  but  is  used  as  a meaningful  excursion  from  the  K-^  case  for  compari- 
son purposes.  An  explicit  representation  using  gaussian  rods  can  be  found  for 
this  case.  Again  let  us  assume  that  there  are  N/A  gaussian  rods  distributed  on 
the  average  over  the  propagation  region.  The  electron  density  is 


n(z.x)  . ^ (2-z,)=^]/a- 


i=l 


where  the  rod  size  and  n^  are  constant.  Let 


(13) 
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^ _ / 4An^A 
° VnL^tt  . 


a = L//2 


N is  distributed  as  in  equation  (4).  The  mean  electron  density  is  found  from 
equations  (13)  and  (14). 


n(z,xT  = TTL2An2j 


The  two-dimensional  correlation  function  of  this  representation  is  exactly 
equation  (12).  The  two-dimensional  power  spectral  don.’ty  is  found  by  trans- 
forming equation  (12). 


PSD  (K  ,K  ) = TiL^An^  ^ K^^-)L2/4 

g X z' 


Again,  neither  the  correlation  function  in  equation  (12)  nor  the  power  spectral 
density  in  equation  (17)  is  a function  of  N/A.  For  this  case  one  MFP  is  de- 


fined as 


(¥)■' 


The  correspondence  between  a MFP  for  this  case  and  a MFP 


for  the  K"2  case  is  evident. 


The  last  correlation  function  corresponds  to  a random  distribution  of  uni- 
form square  striations  with  one  side  perpendicular  to  the  LOS.  The  striations 
are  L on  an  edge.  This  case  was  generated  to  examine  the  effects  of  sharp 
edges  on  the  question  of  determining  when  third  order  or  higher  fluctuation 
statistics  are  important.  The  two-dimensional  power  spectral  density  and  cor- 
relation function  for  this  case  are 


16An^  s 
L2 


in2  \ 2 / sin2  \ 2 / 


I 

L 


R5q(e.p)  = An2  (1 


|/L)(1 


aWBBBWPaE.-l:- I 


AFWL-TR-76-304 


SECTION  III 
PROPAGATION  METHODS 

Three  different  propagation  methods  will  be  discussed.  The  first  two  will 
be  referred  to  as  the  Tatarskii  multi  scatter  (TM)  method  and  the  multiple  phase 
screen  (MPS)  method.  These  two  methods  are  complementary  in  the  sense  that 
there  are  cases  where  only  one  or  the  other  is  practical  and  cases  where  both 
are  applicable.  The  use  of  two  methods  in  regions  where  both  are  applicable 
allows  comparisons  to  check  each  method  respectively  to  insure  accurate  compu- 
tation. The  third  method  called  the  Rytov  approximation  (RA)  method  is  useful 
not  so  much  in  actual  calculation  as  in  the  physical  insight  it  gives  particu- 
larly in  the  onset  of  amplitude  and  phase  flucuations. 

Tatarskii  Multiscatter  Propagation 

The  TM  method  has  been  previously  described  (ref.  7)  and  the  following  has, 
in  part,  been  extracted  from  that  discussion.  Let  the  signal  incident  on  the 
slab  be 


E(z,x)  = (E{z)  + Ep(z,x)  + i Ej(z,x))e^'^^  ^2i) 

where  E(z)  = mean  field 

Ej^(z,x)  = in-phase  quadrature  component 
Ej(z,x)  = out-of-phase  quadrature  component 
z = coordinate  along  LOS 
X = coordinate  perpendicular  to  LOS 
K = field  wave  number 

The  initial  conditions  on  the  field  quantities  are 

EWf  = 1.0  (22a) 


Ep(0,x)  = Ej(0,x)  = 0.0 


(22b) 
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Let  us  now  define  the  complex  autocorrelation  function  and  the  asynmetry 
autocorrelation  functions,  henceforth  to  be  known  as  the  G and  F functions  re- 
spectively. 


G(z,p)  = + E|^(z,x)Ej^(z,x  + p)  + Ejlz,x)Ej(z,x  + p) 


(23) 


F(z,p)  = E(2)2  + Ep(z,x)E(^(z,x  + p) 


- E^{z,x)Ej{z,x  + p)  + 2i  Ep{z,x)Ej{z,x  + p) 


(24) 


since  E^TzTxT  = Ej(z,x)  = 0 and  E|^(z,x)Ej{z,x  + p)  = Ej(z,x)Ej^{z,x  + p).  Tatarskii 
(ref.  2)  derived  the  equation  of  motion  for  G.  The  derivation  for  F follows  in 
a similar  manner  involving  the  same  assumptions.  The  equations  of  motion  are 


dF(z,p) 

dz 


^ - K2F(z,p)[A(z.O)  + A(z.p)] 

dp2 


(25) 


- K2g(z,p)[A(z.O)  - A(z,p)] 


(26) 


where  A(z,p) 


An(z,x)An(z' ,x  + p)  dz' 


(27) 


An  = n - n 
e = electron  charge 
m = electron  mass 
c = speed  of  light 

The  limitations  on  equations  (25)  and  (26)  are  discussed  in  Tatarskii  (ref. 

8)  and  are  summarized  on  the  following  page. 
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X « Jl 


> MFP 


(28a) 


(28b) 


Xa  « 1 


(28c) 


« 1 


(28d) 


aZ  « 1 
S V 


(28e) 


L « X 


(28f) 


where  x = field  wavelength 

a = inner  fluctuation  scale  size 

L = outer  fluctuation  scale  size 

ct  = mean  field  attenuation  coefficient 

= mean  square  scattering  angle 

X = backscattering  attenuation  distance 

These  limitations  do  not  affect  present  work.  Another  implicit  assumption  is 
that  A(z,p)  varies  slowly  with  respect  to  z over  the  largest  scale  size  of  the 
fluctuations.  The  TM  method  is  basically  limited  to  two-dimensional  problems 
except  where  the  fluctuations  are  nearly  isotropic  around  the  propagation  line 
of  sight. 

A propagation  calculation  using  the  TM  method  consists  of  the  following. 
First  the  A(z,p)  function  is  generated.  Next,  equations  (25)  and  (26)  are 
solved  in  the  particular  geometry  desired.  The  F and  G functions  along  with 
the  gaussian  assumption  on  the  quadrature  components  results  in  a complete 
statistical  description  of  the  propagated  signal  structure.  These  statistics 
are  converted  to  time  statistics  by  the  transformation 


= p/v 
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where  7 is  an  average  velocity  calculated  by  considering  the  various  sources  ot 
time  variation  in  the  system.  Finally  sample  signal  structures  are  generated 
from  the  time  dependent  statistics  using  fourier  techniques  described  in  refer- 
ence 9. 

In  addition  to  single  frequency  structures,  the  TM  method  can  calculate  fre- 
quency coherence  bandwidths.  The  coherence  function  is  defined  as 

H(2,p)  = E*  (z,x  + p) 

1^1  1^2  (30) 

H(z,p)  = (z,x)Ep|^^  (z,x  + p)  + Eji^^  {z,x)Ej|^^  {z,x  + p) 


+ i 


( E^j^^'TzTxTEjjJ^TzTxT^  - EjKj  ^ 


(31) 


where  Kj  = first  frequency  wave  number 
K2  = second  frequency  wave  number 

and  H(z,p)  contains  the  mean  field  component  terms.  The  real  part  of  H(z,p)  is 
defined  as  the  frequency  correlation  coefficient.  The  imaginary  part  is  the 
mean  of  the  sine  of  the  phase  difference  excluding  dispersion  between  the  two 
frequencies . 

The  equation  of  motion  for  H(z,p)  is 


dH(z,p) 

dz 


( K]  - K2 ) 
^ 2K1K2 


d2H(z.p) 

dp^ 


/(Ki  + K2)  \ 

41^  \ 2K1K2 


H(z,p) 

(32) 


where  The  form  of  the  equation  is  similar  to  equation  (25)  and 

the  same  solution  technique  can  be  used.  The  initial  condition  is 


H(0,p)  = 1 (33) 

The  coherence  band  width  is  defined  as  the  value  of  |Ki  - K2IC/2  = |vi  - V2I 
where 
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Real  [H(Zg,0)]  = 1/e  = 0.368 

This  coefficient  is  an  estimate  of  the  amount  of  band  spreading  required  to  get 
some  measure  of  frequency  diversity. 

The  TM  method  has  both  significant  advantages  and  disadvantages.  The  basic 
advantages  are  simplicity,  efficiency,  and  a large  range  of  application.  This 
method  can  handle  much  more  severe  conditions  than  the  MPS  or  similar  methods. 
The  main  disadvantage  is  that,  for  realistic  propagation  cases,  it  handles  only 
either  very  weak  scattering  or  Rayleigh  multiscatter  cases.  These  restrictions 
are  a result  of  assuming  gaussian  statistics  for  the  quadrature  components.  It 
is  possible,  in  principle,  to  calculate  the  higher  statistical  moments  using  TM 
techniques  and  do  away  with  the  Rayleigh  limitation,  but,  in  practice,  such 
calculations  are  too  expensive  to  be  practical.  In  addition,  the  gaussian 
assumption  on  the  quadrature  components  is  required  for  efficient  generation  of 
sample  signals  from  the  calculated  statistics.  These  disadvantages  are  not 
very  important  as  the  MPS  method  handles  these  regions  where  the  TM  method  is 
inadequate.  In  addition,  it  will  later  be  apparent  that  the  Rayleigh  scatter 
case  is  the  most  common  case  for  severe  UHF  environments. 

Multiple  Phase  Screen  Propagation 

The  field  satisfies  the  following  scalar  wave  equation. 


{v2  + K2n2(z,x))  E(z,x)  = 0 


where  n^. (z,x)  = n^  + An^. (z,x)  = index  of  refraction 


n 


i 


A . 47re2?P|  ^/2 
\ mc^K^/ 


IT  = mean  electron  density 
e = electron  charge 
m = electron  charge 
c = light  speed 
K = Ztt/X 

X = field  wavelength 


(35) 
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Making  the  usual  parabolic  approximation,  we  have 


[■57--^—-  iKn.  An.(2,x)]  E'(z,x)  = 0 
2Kn.  3x2  1 1 


E’(z,x)  = E(z,x)e'’'^i'^^ 


(36) 

(37) 


If  An^. (z,x)  = 0,  the  field  propagation  is  exactly  soluable. 


00 

E^(z+az,x)  = ^ ^ 


1 4Z  .WttX 

- \T  J e 


2Kn 


ni=-oo 


(38) 


where 


■/: 


■jmiTX 

6^,(2)  = / E‘  (z,x)  e ^ dx 


(39) 


and  T > 5L. 

Equation  (38)  can  be  looked  at  as  a sum  of  plane  waves  of  wave  number  K in- 
dexed by  m traveling  at  an  angle  of  -5^  with  respect  to  the  z direction.  The 

TKn^. 

angles  represented  by  are  always  very  small  for  any  m where  C (z)  is  sig- 

TKn. 

nificant.  This  means  that  the  distance  traveled  by  each  plane  wave  is  very 
close  to  Az  and  hence  each  plane  wave  is  affected  in  approximately  the  same  way 
by  An.(z,x).  If  az  < z where  z,  is  the  distance  a plane  wave  of  wave  number  K 
can  propagate  without  developing  amplitude  fluctuations,  then  the  effects  of 
the  index  of  refraction  fluctuations  on  each  of  the  m waves  is  to  shift  the 
phase  at  z+Az  by 


<(i(x)  = Kn 


An^ (z' ,x)dz' 


(40) 
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The  final  propagated  field  at  z+az  is  approximately 


iil'(x) 


E' (z+Az,x)  = Ej;  (z+AZ,x)e 


Equation  (41)  is  the  basis  of  the  multiple  phase  screen  propagation  algorithm. 
The  striated  layer  of  thickness  z^  is  divided  up  into  layers  of  thickness  az 
< z . The  propagation  of  the  field  through  each  successive  layer  is  done 

“ a 

according  to  equation  (41).  The  determination  of  z will  be  discussed  later  in 
the  discussion  of  the  RA  method.  The  MPS  method  is  made  practical  primarily  by 
the  use  of  the  fast  fourier  transform  algorithm  (FFT)  to  evaluate  equations  (38) 
and  (39)  because  the  FFT  allows  for  a much  more  rapid  transform  than  equation 
(39)  suggests.  The  FFT  produces  a discrete,  evenly  spaced  angular  spectrum, 
which  places  another  limitation  on  az.  It  is  necessary  that 


where 


(m  irN 


is  the  largest  significant  transverse  field  wave  number. 


The  calculation  of  the  integrated  phase  can  be  done  in  two  ways.  The  first 
is  explicit  integration  as  shown  in  equation  (40)  which  requires  that  An^(z,x) 
be  specified.  The  second  way  is  to  determine  (}i(x)  by  Monte  Carlo  techniques 
assuming  that  An^.  is  a gaussian  distributed  variable.  From  equation  (40) 


dz"  An . (z' ,x)An . (z" jX+p] 


Let  z"  = a'  + e 


= W?K2 


J^Z+AZ-Z' 

L. 


dc  R. (e,p) 
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where  R^. (e,p)  = An^.  (z,x)  An^.{2+e,x+p)  is  the  correlation  function  of  the  inoe 
of  refraction  fluctuations.  If  Az  is  greater  than  several  fluctua'-ion  correla- 
tion lengths,  such  that  R^(az,p)  = 0,  then 


4>(x)4>(x+p)  = AZ  n? 


R^.(c,p)de 


(45) 


Since  An^(z,x)  was  assumed  a gaussian  distributed  variable,  (i>(x)  will  be  a 
gaussian  distributed  variable  with  the  above  correlation  function.  Let  (()(x)  be 
represented  as  a fourier  series  on  the  interval  2T. 


m=-oo  (46) 

Since  <t>(x)  is  gaussian  and  zero  mean,  then  so  are  the  b^.  The  variance  of  the 

b„  is  easily  shown  to  be 
m ■' 


—7 — VTT — ; — r i inirx  6 • dp 
<^{x)<(,(x+p)  e — mj 


(47) 


Also  since  4'(x)  is  real,  b^  = b*^.  The  generation  of  a realization  is  now 
rather  easy.  The  fourier  coefficients,  b^,  are  sampled  from  their  known  vari- 
ance. Then  FFT  is  used  to  evaluate  equation  (46)  and  4>(x)  is  then  used  in  equa- 
tion (41).  The  requirement  that  Az  be  several  fluctuation  correlation  lengths 
also  means  that  successive  calculations  of  <t>(x)  will  be  statistically  indepen- 
dent. 

The  MRS  method,  like  the  TM  method,  has  several  advantages  and  limitations. 
The  most  important  advantage  is  that  the  MRS  method  can  handle  the  region  in 
propagation  space  that  includes  weak  scattering  to  near  Rayleigh  scattering 
which  is  the  region  that  the  TM  method  cannot  handle.  A second  advantage  is 
that  when  the  FFT  is  used,  the  MRS  is  relatively  fast  and  inexpensive  so  many 
cases  can  be  run.  The  limitations  on  the  MRS  primarily  stem  from  numerical 
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details.  First,  the  method  is  practical  only  for  two-dimensional  problems  due 
to  computer  considerations.  This  is  not  a very  serious  problem  as  propagation 
before  the  Rayleigh  limit  is  reached  is  primarily  two-dimensional.  The  most 
serious  limitation  on  the  MRS  method  derives  from  the  practical  limitation  on 
the  smallness  of  Ax,  the  transverse  mesh  spacing.  As  the  propagation  environ- 
ment degrades  the  distance  over  which  significant  phase  or  amplitude  changes 
occur  becomes  smaller.  In  practice,  it  is  required  that 


AX  < : [L/( K^An?  Zg)] 


This  restriction  limits  the  severity  of  effects  that  can  be  handled  by  the  MRS 
method.  Fortunately,  calculations  show  that  UHF  signals  become  Rayleigh  before 
AX  becomes  so  small  as  to  make  MRS  calculations  impractical.  The  TM  method  can 
be  used  in  the  Rayleigh  regime  after  the  MRS  fails. 

Rytov  Approximation  Method 

The  field  must  satisfy  the  scalar  wave  equation  in  equation  (35). 


ij;(z,x)  = log  (E(z,x)) 


v2(()  + + K^n?(z,x)  = 0 


>1^(2. x)  = <^q(z,x)  + i|;i(z,x) 


n^.(z,x)  = h^.  + An^.  (z,x) 
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where 


^^(z,x)  and  n.  satisfy  equation  (50)  for  the  unperturbed  case. 


+ (vi1)q)2  + n?  = 0 


(53) 


If 


lAn.(z,x)l  « 1 


(54) 

(55) 


Then 


+ 2viJ|q  • An^(z,x)  “ 0 

If  = X ■*■  ■'<(>»  after  much  manipulation  (ref.  8)  it  can  be  shown 


(56) 


fe  / 

-CO 


(z,e,Kp)e^^  ’^p^ 


(57) 


dz  cos2[^J  (z^+Zg-z)^  der. 


,.(z,e,K)e  p 


(58) 


e"^ 

:.(z,e,K)  =1  An.(z,x)An^(z+e,x+p) 

>-00 


-iK  p j 
p dp 


(59) 
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The  z dependence  in  C^(z,e,Kp)  represents  changes  in  the  index  of  refraction 
statistics  along  the  propagation  line  of  sight.  This  method  has  all  the  basic 
restrictions  of  the  TM  method.  In  addition,  it  is  required  that  s 0-3.  The 
basic  value  of  this  method  is  that  it  gives  insight  into  the  onset  and  devel- 
opment of  amplitude  fluctuations,  and  phase  fluctuations  due  to  propagation. 

Examination  of  equations  (57)  and  (58)  provides  a prescription  for  deter- 
mining when  amplitude  fluctuations  occur  and  their  dependence  on  the  environ- 
ment properties.  First  if 


L^K 


C.(z,e,K  ) » C,(0,e,K  ) 

1 P 1 P 


(60a) 

(60b) 


then 


~ /2ne^^  An^ 
\mc2  / k2 


(61) 


^2  ; / 2neiV  jt 

\mc2  / 


(62) 


For  this  case,  the  amplitude  and  phase  fluctuations  occur  simultaneously.  This 
is  typical  of  VHF  satellite  links.  If  condition  (60a)  is  not  true  then  the  re- 
sults depend  on  the  particular  power  spectrum  assumed.  When 


L^K 

z+2g 


» 1 , 


X2  . 1 

^ T2 


[(Zg-<-Z^)3-2^3] 

K** 


&)’/'  '■ f 


dK 


(63) 


For  a gaussian  spectrum,  equation  (63)  is  exactly  calculable. 
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2 = ^ 


K'*L3 


Vtnc^  / 


(6‘+, 


For  the  K"2  spectrum,  things  get  more  complicated  since 


/ 


C(e,K  )de  ; 

P (1/L2  + K2)3/2 


(65) 


until  the  inner  scale  is  reached.  This  means  that  the  integrand 
in  equation  (63)  increases  until  the  spectrum  inner  scale  is  reached  or  until 

K 


After  some  manipulation 


= C 


\mc2  / 


K3L 


^s'*‘^g 


(66) 


where  Cj  is  a constant  that  was  found  by  direct  calculation  to  be  approximately 
equal  to  0.5. 

Let  us  summarize.  For  the  K~2  spectrum  we  have 

p 


= (?Miy  !si,  ^ , 


\mc2  / K2 


s 9 


(67a) 


= /?Mif  0.5 > 

\mc2  / LK^  ^s  ^g 


(67b) 


For  the  gaussian  spectrum  we  have 
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The  initia]  phase  fluctuations  go  as 


(69a) 


(69b) 


for  either  spectrum. 

In  general,  equations  (67b),  (68b),  and  (69b)  are  usually  appropriate.  By  taking 
the  ratio  of  equation  (68b)  to  equation  (67b)  we  get 


For  — — >>  1,  the  spectrum  results  in  a quicker  onset  of  amplitude  fluctu- 
ations.  Taking  the  ratio  of  equations  (67b)  and  (69b)  we  get 


which  means  that  the  phase  fluctuations  onset  before  the  amplitude  fluctuations. 
This  is  observed  in  ambient  environments  and  in  the  calculations  to  follow. 

These  equations  show  the  advantages  of  maximizing  K. 
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SECTION  IV 
RESULTS 


The  calculations  to  be  discussed  all  assumed  a nominal  UHF  frequency  of  300 
Mhz.  The  initial  field  strength  and  field  amplitude  was  unity. 

Figures  2 through  7 show  several  quadrature  component  correlation  functions 
for  the  and  intermediate  spectrums.  Several  of  these  figures  compare  the 
TM  and  MRS  results  for  identical  scenarios  and  show  good  agreement.  The  agree- 
ment for  the  strong  scatter  cases  is  within  the  statistical  noise  of  the  MPS 
calculation.  When  the  signal  structure  is  Rayleigh,  these  correlation  func- 
tions completely  define  the  signal  statistics  and  sample  signals  can  be  gener- 
ated using  methods  discussed  in  reference  7.  Even  when  the  statistics  are  not 
Rayleigh,  these  correlation  functions  are  very  useful  because  their  width  is  a 
measure  of  the  signal  structure  correlation  distance  when  the  mean  square  phase 
is  greater  than  about  one.  Let  us  define  the  correlation  distance  (n  ) as  the 
e"^  point  of  the  correlation  function  curve.  Figure  8 shows  as  a function 
of  the  environment  parameters  for  the  and  the  intermediate  spectrums.  In 
figure  8 a distinction  is  made  between  an  outer  scale  size  parallel  to  the 
propagation  line  of  sight  (Ln)  and  a perpendicular  outer  scale  size  (L4_  = 

1.25  X 10^).  This  separation  is  necessary  to  parameterize  the  behavior  of  j!.^. 
For  layers  with  uniform  statistical  properties  can  be  adequately  represented 
by  a power  law  function. 

) 

^0  “Lj.(Lu  An2  Zg/K2)-fi  ^72) 


For  nonuniform  layers  the  can  be  expressed  as  a power  law  function  of  an 
appropriately  defined  integral.  Let  us  define  for  the  and  the  intermediate 
spectrum  an  integral  called  the  scale  size  factor. 


SSF 


An2  dz 

K2Lj.'-"2  3 


(73) 
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Then 


dz 

K2Lj_1.941 


,,  -0.616 
I = 2.4  X IQi**  SSF„.2 
0 K ^ 


(74) 


(75) 


1 


I = 1.24  X 1012  sSF;0.515 
0 Inr. 


(76) 


The  fit  for  the  K"2  spectrum  assumes  that  the  environment  inner  scale  size 
is  zero.  The  curve  set  by  the  X-s  assumed  a 10  meter  inner  scale  and  hence  the 
curve  breaks  upward  as  = 10  meters  is  approached.  The  slope  of  the  "X" 
curve  approaches  minus  one  half  asymptotically  as  the  scale  size  factor  in- 
creases. 

If  <v>  is  the  velocity  of  the  line  of  sight  across  the  striations,  then  the 

correlation  time  x can  be  defined. 

0 


To  = y<v> 


(77) 


Thus  equations  (73)  through  (76)  with  <v>  define  for  the  two  discussed  spectra 
the  basic  time  dependence  of  the  signal  parameters  for  mean  square  angles 
greater  than  one.  Examination  of  figures  4 through  7 also  shows  that  the  cor- 
relation functions  are  approximately  gaussian  in  shape  so  only  one  parameter  is 
necessary  to  define  the  second  order  signal  statistics.  Thus  for  Rayleigh  sta- 
tistics, the  statistical  description  is  complete  with  one  parameter.  This  is 
an  important  result  in  that  it  implies  that  in  Rayleigh  environments  system  per- 
formance is  only  a function  of  x^  and  the  signal-to-noise  ratio. 

For  the  case  where  the  signal  disturbances  are  mainly  phase  fluctuation  with 
little  or  no  amplitude  fluctuation,  can  be  related  to  the  phase  correlation  ,, 

function.  Let  J 

}' 
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(P(x)<p(x+p)  - <p^  F(p) 


(7S' 


If  <t>(x)  is  gaussian  distributed,  as  postulated  for  thick  layers,  then  it  can  be 
shown  from  equation  (2F)  that 


E*(x)E(x+p)  = E2  e”'^^  (1-F(p)) 


(79) 


Then 


F(«„)  = 1 - L 

a2 


(80) 


when  (j|2  >>  1.  Equation  (78),  which  defines  F(p),  and  equation  (80)  demonstrate 


a clear  link  between  a and  the  phase  fluctuations.  If  we  assume 

0 


F(p)  - 1 - P^/Pq^  + — 


E*(x)E(x+p)  ==  1 - p2/z^2  + 


(81) 


(82) 


then 


Po^  “ 


(83) 


Equation  (83)  demonstrates  the  basic  relationship  between  p^,  the  phase  correla- 
tion distance,  and  Equation  (83)  also  shows  that  the  phase  correlation 
distance  is  proportional  to  the  root-mean-square  phase  fluctuation,  a fact  ig- 
nored in  many  works.  Let  us  estimate  the  rate  of  change  of  the  phase  since 
this  is  the  important  quantity. 


1/2 


d 

dx 


—1/? 
d)2  5, 


(84) 


27 


AFWL-TR-76-304 


This  shows  that  the  mean  square  phase  fluctuation  which  is  often  used  to  esti- 
mate possible  system  degradation  in  reality  contains  little  useful  information 
without  accompanying  knowledge  of  L or 

Figures  9 through  13  show  the  quadrature  component  distribution  functions 
for  the  K“2  spectrum  as  Ana  increases.  The  solid  curves  in  the  last  two  compo- 
nent distribution  functions  are  gaussian  curves  which  show  the  distribution 
limit  as  Ana  becomes  large.  These  plots  show  that  the  amplitude  fluctuations 
become  significant  with  only  a few  radians  of  root-mean-square  phase  fluctua- 
tion. Indeed,  as  the  component  correlation  functions  become  equal  the  ampli- 
tude fluctuations  result  in  a near  unity  scintillation  index.  The  statistics 
become  essentially  gaussian  when  Ana  ==  10^°. 

Figures  14  through  17  and  18  through  23  show  the  component  distributions 
for  the  intermediate  and  gaussian  spectrum  respectively.  The  results  for  the 
intermediate  spectrum  is  little  different  from  the  gaussian  so  only  the  latter 
will  be  discussed.  As  with  the  spectrum,  the  component  distributions  be- 
come equal  after  a few  radians  of  root-mean-square  fluctuation.  The  amplitude 
fluctuations,  however,  develop  more  slowly  with  a different  character.  The 
population  of  deep  fades  is  less  than  with  the  K'^  spectrum  up  to  An2  = 
and  then  exceeds  the  spectrum  before  settling  down  to  a Rayleigh  distribu- 
tion at  An2  = 10^2.  The  population  of  a signal  enhancements  follows  no  simple 
behavior.  The  sides  of  the  component  distributions  for  the  gaussian  spectrum 
are  steep,  dipping  below  the  gaussian  distribution  but  going  above  for  very 
large  amplitudes.  This  behavior  is  simple  to  understand.  For  the  gaussian 
spectrum  there  is  essentially  a single  size  of  fluctuations  at  about  K’l  = 1.25 
km.  The  receiver  at  z^  = 200  km  is  in  the  near  zone  until  Anz  becomes  very 
large.  The  initial  propagation  effects  manifest  themselves  as  large  phase 
fluctuations  with  only  small  amplitude  fluctuations.  Eventually,  amplitude 
fluctuations  develop  as  Anz  increases.  The  K"^  spectrum  has  the  same  large 
size  as  the  gaussian  but  also  has  a significant  population  of  smaller  sizes 
whose  zone  the  receiver  is  in  or  near.  This  gives  rise  to  amplitude  fluctua- 
tions at  much  smaller  fluctuation  power  than  does  the  gaussian  spectra. 

Additional  insight  is  gained  by  examining  figures  24  through  27  which  are 
selected  amplitude  distributions  for  the  three  spectrums.  These  plots  illus- 
trate the  observations  made  from  the  component  distributions.  The  population  of 
fades  of  a given  depth  is  particularly  apparent.  Initially,  the  K"^  spectrum 


28 


AFWL-TR-76-304 


gives  the  largest  population  of  deep  fades  but  the  gaussian  spectrum  soon 
catches  up.  An  important  point  in  the  K"^  amplitudes  is  that  the  I'opulation  of 
deep  fades  is  never  greater  than  the  Rayleigh  limit.  This  means  L,,at  for  sys- 
tems that  are  primarily  sensitive  to  fades,  the  Rayleigh  limit  is  the  worst 
case  for  a K"^  spectrum.  For  the  gaussian  spectrum  the  situation  is  quite 
different.  Let  us  assume  that  the  fades  are  slow  and  that  a system  has  10  dB 
link  margin.  The  performance  will  be  a function  of  the  population  of  fades 
with  amplitudes  less  than  0.316.  It's  clear  from  figure  27  that  for  An2  > 3.2 
X 10^°,  the  gaussian  spectrum  will  result  in  a poorer  performance  than  the 
Rayleigh  limit. 

The  difference  between  the  and  the  gaussian  spectrums  can  be  seen  in 
figure  28  which  contains  the  six-bit  character  error  rate  curve  for  the  Tacti- 
cal Transmission  System  (TATS)  wideband  modem.  This  system  is  a 8-array  Fre- 
quency shift  key  UHF  modem  that  was  studied  in  detail  in  a previous  work  (ref. 
9).  The  same  model  was  used  here.  The  solid  curve  is  from  the  K"^  spectrum. 

For  this  spectrum  the  statistics  are  Rayleigh  for  t < 0.5  seconds.  The  re- 
mainder of  the  curve  reflects  the  system  performance  as  the  Rayleigh  limit  is 
approached  and  is,  in  general,  a function  of  more  than  t . The  important 
point  to  note  is  that  in  the  approach  to  Rayleigh  statistics,  the  performance 
was  always  better  than  that  associated  with  the  Rayleigh  limit  confirming  the 
earlier  discussion  about  the  amplitude  distribution.  For  the  gaussian 
case,  the  Rayleigh  limit  is  not  reached  until  t < 3 x 10'^  seconds  at  which 
point  the  two  curves  are  nearly  convergent.  Sine  the  curve  represents  the 

p 

Rayleigh  limit,  it's  clear  that  the  gaussian  spectrum  gives  performance  which 
is  worse  than  Rayleigh  as  was  suggested  earlier. 

The  basic  issue  for  power  spectra  in  UHF  propagation  is  whether  real  spec- 
trums are  most  often  multisize  such  as  the  K"^  or  single  size  like  the  gaussian. 
If  is  the  rule,  then  Rayleigh  statistics  with  given  pose  a reasonable 
worst  case  threat  to  UHF  systems.  At  least  ninety  percent  of  ambient  scintil- 
lation measurements  have  been  K'^,  so  reasonable  worst  case  analysis  is  pos- 
sible. For  other  disturbed  environments,  the  correct  spectrum  or  spectrums  are 
not  yet  well  defined.  This  study  suggests  that  emphasis  should  be  put  on  re- 
solving the  power  spectrum  unknowns. 
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The  last  question  to  be  examined  is  the  affect  of  third  order  or  higher 
striation  statistics  on  the  propagation  results.  Figure  29  contains  the  ampli- 
tude distributions.  Si,  and  £ as  a function  of  the  mean-free-path  thickness  pf 

0 

the  scattering  layer  for  a gaussian  power  spectral  density.  The  infinite  mean- 
free-path  case  represents  gaussian  distributed  fluctuations.  Despite  the  sta- 
tistical jitter  in  the  caluclations  the  thinest  case  is  clearly  different  from 
the  thickest  case.  It  appears  that  five  to  ten  mean  free  paths  are  needed  to 
get  reasonable  insensitivity  to  environment  statistics  above  second  order. 

Figures  30  and  31  show  the  same  propagation  information  as  figure  29,  only 
the  power  spectral  density  is  now  with  the  parameters  as  shown  in  the 
figures.  These  figures  imply  that  three  or  more  mean  free  paths  are  suffi- 
cient to  become  insensitive  to  the  higher  order. 

Figures  32  and  33  repeat  figures  30  and  31  with  the  inner  scale  essentially 
zero  and  including  gaussian  distributed  fluctuations.  Figure  34  repeats  the 
calculations  represented  in  figure  33  with  Anz  reduced  to  3.33  x 10®  and  with 
the  phase  plotted  instead  of  the  amplitude.  Again,  as  above,  only  a few  mean 
free  paths  are  necessary  to  randomize  the  phase  and  become  dependent  only  on 
the  first  and  second  order  fluctuation  statistics. 

Figure  35  contains  results  assuming  square  top  hat  striations.  The  sharp 
edges  scatter  so  effectively  that  three  mean  free  paths  agree  closely  with 
assuming  gaussian  fluctuations. 

It  appears  clear  that,  for  most  applicators,  the  second  order  statistics 
or,  equivalently,  the  fluctuation  correlation  function  is  a necessary  and  suf- 
ficient environment  description  for  propagation  purposes.  The  higher  order 
statistics  can  be  supplied  by  any  reasonable  assumption  such  as  gaussian  dis- 
tributed fluctuations. 
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SECTION  V 
CONCLUSIONS 


This  paper  has  examined  the  propagation  of  electromagnetic  signals  through 
a region  with  fluctuations  in  the  index  of  refraction.  The  region  was  assumed 
to  be  a thick  multiscatter  layer  characterized  by  one  of  four  types  of  power 
spectral  densities.  These  power  spectrums  were  called  the  K"^,  intermediate, 
gaussian,  and  square  top  hat  power  spectrums.  The  calculation  results  indi- 
cated that  the  intermediate  and  the  gaussian  cases  were  sufficiently  similar 
that  only  one  of  them,  the  gaussian,  was  discussed  in  detail. 

The  spectrum  is  by  far  the  most  common  measured  spectrum  in  disturbed 
natural  environments.  The  outer  scale  size  associated  with  the  largest  mea- 
sured fluctuation  power,  Anz/p  = 0.8,  is  L = 1.25  km.  This  spectrum  is  con- 
sidered by  some  to  be  a candidate  for  nuclear  disturbed  environments.  The 
gaussian,  intermediate,  and  square  striation  spectrum  have  also  been  used  for 
disturbed  environment.  These  spectrums  vary  from  the  K~^  spectrum  in  that  they 
are  characterized  by  a single  size.  Examination  of  turbulence  spectrums  re- 
veals that  the  K"^  and  the  gaussian  spectrum  bracket  the  reasonable  range  of 
p 

spectrums.  The  basic  conclusions  concerning  propagation  through  these  spec- 
trums are: 


1.  The  onset  of  phase  and  amplitude  fluctuations  is  determined  by  setting 
^ ^ = 0.1  in  the  following  equations  and  assuming  constant  statis- 

tical properties. 

I 2k 

If  < 1.  then 

s g 


X 
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(67a) 
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L^K 
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> 1 , then 
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(67b) 


(gaussian  spectrum) 
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2.  The  Rayleigh  limit  is  a reasonable  worst  case  for  the  K“2  spectrum  for 
UHF  signals.  This  means  that  the  reasonable  worst  case  performance  of 
UHF  systems  can  be  parameterized  with  only  two  parameters,  and  the 
signal  to  noise  ratio.  The  Rayleigh  limit  is  not  a worst  case  for  the 
gaussian  specturm.  For  this  latter  case,  there  is  a focusing  region 
where  the  population  of  deep  fades  exceed  the  Rayleigh  case.  This  case 
does  not  appear  to  be  significant  because  of  a low  probability  of  occur- 
rence. 


3.  The  distance  {i^)  over  which  the  phase  or  amplitude  can  change  signifi- 
cantly is  specified  by  the  following  equations.  For  the  K"^  spectrum, 
we  get 


SSF 

% 


L|iAn2  dz 

K2l, 1.623 


i = 2.4  X IQi*^  SSF'O-616 
0 

For  the  intermediate  spectrum,  we  get 

^ LiiAn^  dz 

K2l, 1.941 


(73) 

(75) 


(74) 
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Z = 1.24  X 1012  ssf-O-515 
o 

4.  For  both  spectrums,  5 to  10  mean  free  paths  are  sufficient  to  become 
sensitive  only  to  the  fluctuation  correlation  function.  This  conclu 
sion  is  very  important  as  it  implies  that  for  most  cases  of  interest 
only  a limited  amount  of  environment  information  is  required. 


(76) 
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Figure  1.  Propagation  Geometry. 
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Figure  2.  Quadrature  Component  Correlation  Functions  for  An^  = 1.23  x 10^. 
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Figure  3.  Quadrature  Component  Correlation  Functions  for  An^  = IQS. 
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Figure  7.  Quadrature  Component  Correlation  Functions  for  An 


Figure  8.  Field  Correlation  Distance  (t  ) 
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Figure  10.  Quadrature  Component  Distributions  for  K"^  Spectrum  and  An^  = 10®. 


Figure  12.  Quadrature  Component  Distributions  for  K"^  Specturm  and  An^  = 1.23  x 10^. 
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Figure  15.  quadrature  Component  Distributions  for  Intermediate  Spectrum  and  An^  = 1.23  x 10^. 
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Figure  16.  Quadrature  Component  Distributions  for  Intermediate  Spectrum  and  An 
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Figure  17.  Quadrature  Component  Distributions  for  Intermediate  Spectrum  and  An’ 
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Figure  18.  Quadrature  Component  Distributions  for  Gaussian  Spectrum  and  An 
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Figure  22.  Quadrature  Component  Distributions  for  Gaussian  Spectrum  and  An 


Figure  24.  Amplitude  Distributions  for  the  K“2  Spectrum  With  L = 1.25  Km 


Figure  26.  Amplitude  Distributions  for  the  Intermediate  Spectrum. 
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Figure  27.  Amplitude  Distributions  for  the  Gaussian  Spectrum. 
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Figure  30.  /'mplitude  Distributions 
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Figure  33.  Amplitude  Distributions 
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Figure  35.  Square  Top  Hat  Amplitude  Distributions 
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